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Abstract 


Exploring the properties of exoplanets near or inside the radius valley provides insight on the transition from the 
rocky super-Earths to the larger, hydrogen-rich atmosphere mini-Neptunes. Here, we report the discovery of TOI- 
1452b, a transiting super-Earth (R, = 1.67 + 0.07 RK) in an 11.1 day temperate orbit (T,, = 326 + 7 K) around the 
primary member (H = 10.0, Tere = 3185 + 50 K) of a nearby visual-binary M dwarf. The transits were first detected 
by the Transiting Exoplanet Survey Satellite, then successfully isolated between the two 3/’2 companions with 
ground-based photometry from the Observatoire du Mont-Mégantic and MuSCAT3. The planetary nature of TOI- 
1452b was established through high-precision velocimetry with the near-infrared SPIRou spectropolarimeter as 
part of the ongoing SPIRou Legacy Survey. The measured planetary mass (4.8 + 1.3 Mq) and inferred bulk density 
(G652 g cm °) is suggestive of a rocky core surrounded by a volatile-rich envelope. More quantitatively, the mass 
and radius of TOI-1452b, combined with the stellar abundance of refractory elements (Fe, Mg, and Si) measured 
by SPIRou, is consistent with a core-mass fraction of 18% + 6% and a water-mass fraction of 22'7}%. The water 
world candidate TOI-1452b is a prime target for future atmospheric characterization with JWST, featuring a 
transmission spectroscopy metric similar to other well-known temperate small planets such as LHS 1140b and K2- 
18 b. The system is located near Webb’s northern continuous viewing zone, implying that is can be followed at 


almost any moment of the year. 


Unified Astronomy Thesaurus concepts: Exoplanets (498); Super Earths (1655); M dwarf stars (982) 


Supporting material: machine-readable table 


1. Introduction 


Over the past decade, it has become increasingly clear that 
the typical extrasolar planetary system is quite different from 
our solar system. Exoplanets are usually found in a much more 
compact orbital configuration (Howard et al. 2010), and the 
majority of systems have at least one planet with a size 
intermediate between the Earth and Neptune (Howard et al. 
2012; Fressin et al. 2013). Population studies based on the 
Kepler sample have shown that the occurrence rate distribution 
of close-in (P < 100 days) exoplanets displays a valley/gap 
near 1.5—2.0R, (Fulton et al. 2017; Fulton & Petigura 2018; 
Mayo et al. 2018; Hardegree-Ullman et al. 2020). This radius 
valley most likely separates scaled-up, rocky versions of the 
Earth (super-Earths) and hydrogen-rich planets reminiscent of 
Neptune, but smaller (mini-Neptunes). This transition is known 
to be period-dependent (Van Eylen et al. 2018; Martinez et al. 
2019) and to vary with host star properties such as metallicity 
(Owen and Murray-Clay 2018; Petigura et al. 2018), mass 
(McDonald et al. 2019; Cloutier & Menou 2020), and age 
(Berger et al. 2020; David et al. 2021). The existence of a 
radius valley was rapidly attributed to total or partial 
photoevaporation of the atmosphere by highly energetic 
photons during the first 100 Myr, when the host star is more 
active (Owen & Wu 2013; Lopez & Fortney 2014; Owen & 
Wu 2017; Lopez & Rice 2018; Wu 2019). However, another 
atmospheric erosion mechanism is plausible, involving mass 
loss caused by the release of energy from the planet core, 
accumulated during formation and slowly cooling down over 
Gyr timescales (Ginzburg et al. 2018; Gupta & Schlicht- 
ing 2019, 2020). More recently, Lee & Connors (2021) have 
shown that the radius valley can be sculpted as a feature of 
formation, involving gas-poor accretion and supporting the 
hypothesis of a primordial bimodal distribution, rather than the 
result of subsequent atmospheric erosion. In order to identify 
which mechanism dominates, Rogers et al. (2021) predict that 
the number of well-characterized small exoplanets must reach 
25000. Such characterization requires the precise knowledge 
of the planetary radii (<5% uncertainty) and, if possible, the 
planet mass. The combination of the two measurements leads 
to the mean density of the objects, a way to determine whether 


their internal structure is compatible with a rocky, gaseous, or 
intermediate bulk composition. 

Identifying new small planets transiting nearby bright stars is 
the primary objective of the ongoing NASA _ Transiting 
Exoplanet Survey Satellite (TESS) mission (Ricker et al. 
2015). In operation since 2018, TESS has observed 85% of the 
celestial sphere, staring for at least ~27 days at over 50 sectors 
covered so far (24° x 96° per sector). The TESS survey has 
already unveiled more than 5000 candidate exoplanets of 
which more than 200 have been confirmed as new transiting 
planetary systems, including small planets around M-dwarf 
hosts (e.g., TOIJ-270, Giinther et al. 2019; LP 791-18, 
Crossfield et al. 2019; L98-59, Cloutier et al. 2019; 
LTT 1445 A, Winters et al. 2019; LTT 3780, Cloutier et al. 
2020a; TOI-1235, Cloutier et al. 2020b; TOI-700, Gilbert et al. 
2020; TOI-1266, Demory et al. 2020; LP 714-47, Dreizler et al. 
2020; TOI-776, Luque et al. 2021). 

M dwarfs represent prime targets not only for TESS, but in 
exoplanetary science in general. They are the most abundant 
stars in the solar neighborhood (Reylé et al. 2021) and host on 
average 2.5+0.2 planets per M dwarf with radii 14 Ro 
(Dressing & Charbonneau 2015). Their smaller size compared 
to Sun-like stars facilitate the detection and characterization of 
new exoplanets by producing deeper transits for planets of a 
given size. The larger planet-to-star mass ratio amplifies the 
planetary radial velocity (RV) signal, allowing easier mass 
determination. Lastly, their lower luminosity results in a closer- 
in habitable zone (HZ), with orbital periods typically of one or 
two weeks adequately sampled by TESS. 

The James Webb Space Telescope (JWST) is poised to 
revolutionize the field of exoplanet atmospheres (Bean et al. 
2018) by offering a collecting area more than 6 times larger 
than the Hubble Space Telescope (HST) and spectral coverage 
from the visible to the mid-infrared (0.6-28 wm). JWST will 
allow simultaneous identification of many chemical species 
with large absorption cross section in the infrared (e.g., H2O, 
CHy,, CO, COz, NH3), as well as probe the atmosphere of 
terrestrial planets with unprecedented sensitivity. One key 
objective of TESS is to discover the best transiting exoplanets 
amenable for atmospheric characterization with JWST (Kemp- 
ton et al. 2018). 
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TESS Target Pixel File (21”/pixel) 
4 A OU —. 


Table 1 
TESS Observations of TOI-1452 
Sector Camera CCD UT Start Date UT End Date 
14 3 2 2019-07-18 2019-08-14 
15 3 2 2019-08-15 2019-09-10 
16 2 1 2019-09-12 2019-10-06 
17 4 2 2019-10-08 2019-11-02 
19 4 1 2019-11-28 2019-12-23 
20 4 1 2019-12-24 2019-01-20 
21 4 1 2020-01-21 2020-02-18 
22 4 4 2020-02-19 2020-03-17 
23 4 4 2020-03-19 2020-04-15 
24 3 4 2020-04-16 2020-05-12 
25 3 3 2020-05-14 2020-06-08 
26 3 3 2020-06-09 2020-07-04 
40 3 2 2021-06-25 2021-07-23 
41 3 2 2021-07-24 2021-08-20 
47 4 1 2021-12-31 2022-01-27 


Here, we report the discovery of a new small exoplanet around 
the nearby M dwarf TOI-1452. The planet was first detected by 
TESS, then characterized via follow-up efforts including RV 
monitoring with the SPIRou spectropolarimeter. The complete set 
of observations is described in Section 2. The host star properties 
and physical parameters are derived in Section 3. Our data analysis 
and results are presented in Section 4. The implications of this 
discovery and prospects for follow-up characterization are discussed 
in Section 5, followed by concluding remarks in Section 6. 


2. Observations 
2.1. TESS Photometry 


TOI-1452 (TIC 420112589) was observed by TESS in sectors 14 
through 26 (except 18), thus almost continuously from 2019 July 18 
to 2020 July 4, in sectors 40-41 in 2021 from June 25 to August 20, 
and finally in sector 47 from 2021 December 31 to 2022 January 27 
(details in Table 1). TOJ-1452 was sampled at the TESS 2 minutes 
“short” cadence, as the star is part of the cool dwarf list (Muirhead 
et al. 2018), a specially curated list of high-priority late-K and M 
dwarfs added to the TESS Input Catalog (TIC; Stassun et al. 
2018b, 2019). We used the publicly available*® per sector light 
curves produced by the TESS Science Processing Operations 
Center (SPOC; Jenkins et al. 2016) at NASA Ames, more 
specifically their Presearch Data Conditioning Simple Aperture 
Photometry (PDCSAP; Smith et al. 2012; Stumpe et al. 
2012, 2014). The PDCSAP light curves are corrected for both 
instrumental systematic trends seen across stars in the same 
sector/camera/CCD and for flux contamination from nearby 
stars located within a few TESS pixels (21”). Flux dilution 
reduces the observable transit depth, resulting in an under- 
estimation of the planetary radius if not accounted for. This was 
particularly important for TOI-1452 because a companion star 
(TIC 420112587; see Section 3.2) is separated by only 3”2 and 
has a similar magnitude in the TESS band (AT = 0.204). A new 
background correction was implemented for the TESS extended 
mission (starting with sector 27). We followed the procedure 
outlined in the TESS Data Release 38 notes** to correct our 
PDCSAP fluxes from the primary mission (sectors 14-26), 


43 Mikulski Archive for Space Telescopes (MAST): archive.stsci. 
edu/tess/. 
archive.stsci.edu/tess/tess_drn.html 
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Figure 1. Upper panel: TESS target pixel file of TOI-1452 from sector 14 
(11 x 11 pixels subregion). The optimal aperture to extract the PDCSAP 
photometry is delimited by a black line. Nearby Gaia EDR3 sources are 
represented with blue circles. TOI-1452 and its 3/2 companion (TIC 
420112587) are shown with a red and a cyan circle, respectively. Lower 
panel: Same region of the sky observed with OMM-PESTO 1.6 m on 2020 
February 22 (see Section 2.2). The visual binary was partially resolved. 


adjusting the baseline level and reducing the inferred transit 
depth by ~1.7%. This ensures that the primary and extended 
mission data produce the same estimate of the planetary radius. 
Figure | shows an 11 x 11 pixel subregion around TOI-1452 
from TESS sector 14 and the same region of the sky observed 
from the ground. This illustrates how TESS alone cannot resolve 
the source of a transit between TOI-1452 and the nearby 
companion TIC 420112587. The normalized PDCSAP light 
curve of TOI-1452 from sectors 14 and 21 is presented in 
Figure 2, while the remaining sectors are shown in Figure A1. 

A search of the sectors 14-16 with an adaptive, wavelet- 
based matched filter (Jenkins 2002; Jenkins et al. 2010, 2020) 
first identified transit signatures for TOI-1452. The data 
validation reports (DVR; Twicken et al. 2018; Li et al. 2019) 
fitted a limb-darkened transit model with a signal-to-noise ratio 
(SNR) of 8.0, a period of 11.06409 days, and an average 
uncontaminated transit depth of 3.778 parts per thousand (ppt), 
corresponding to a preliminary planetary radius of 
1.83 + 0.30 Rz. This led to the announcement of the planet 
candidate TOI-1452.01 (Guerrero et al. 2021) by the TESS 
Science Office on 2019 October 26. Simultaneously, the TESS 
mission announced the candidate TOI-1760.01 around the 
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Figure 2. Left panels: Normalized PDCSAP light curve of TOI-1452 from sectors 14 and 21, featuring transits (blue data points), a ~5% stellar flare event (zoomed in 
subpanel), and outliers (red data points) either rejected by sigma clipping (3.50 clip) or manually (sector 21). A quasiperiodic Gaussian process model is depicted with 
the green curve (details in Section 4.2). The remaining sectors are presented in Figure Al. Right panel: TESS phase-folded corrected transits (32) from sectors 14-26, 
40-41, and 47. Binned photometry (8 minute phase bin) is represented with black points. The blue curve shows the best-fit transit model (described in Section 4.3), 
with the 68% confidence interval envelope in light blue. The residuals of this fit are shown below. 


companion star TIC 420112587 sharing the same ephemeris as 
TOI-1452.01. Ground-based photometry was able to isolate the 
transit signal, originating from TOI-1452 (see Section 4.1). The 
latest available DVR from sectors 14—41 includes 30 transits 
and reports a period of 11.06196 days along with a radius 
R,= 1.60 + 0.42 R. Our complete reanalysis of the TESS 
light curve presented in Section 4 has resulted in a more precise 
planetary radius, in agreement with previous estimates. 
Figure 2 shows the phase-folded 32 transits from currently 
available sectors. We note that TESS is expected to continue 
monitoring TOI-1452 during 2022. 


2.2. OMM-PESTO Transit Monitoring 


Due to the coarse image sampling of TESS (21” per pixel), 
the origin of a transit signal may be ambiguous when several 
stars are located inside the aperture (e.g., Figure 1). For this 
reason, TESS planet candidates are prone to false positives, 
occasionally attributed to a nearby eclipsing binary (NEB) 
contaminating the light curve (Sullivan et al. 2015). Ground- 
based follow-up with arc-second angular resolution is therefore 
necessary to validate on-target transit and reject the NEB 
scenario. For TOJ-1452, a particular challenge was to 
determine the signal’s provenance between the target and its 
3/2 companion. 

Two transit events of TOI-1452.01 were observed using the 
PESTO camera installed on the 1.6m telescope of the 
Observatoire du Mont-Mégantic (OMM), Québec, Canada. 
PESTO features a 1024 x 1024 pixel EMCCD detector with a 
pixel scale of 0466, providing a field of view (FOV) of 
7195 x 7195. We scheduled the two observing sequences with 
the TESS Transit Finder (TTF), a customized version of 
the Tapir software package (Jensen 2013), and have used 
AstrolImaged (AIJ; Collins et al. 2017) to perform image 


calibrations, including bias subtractions and flat field correc- 
tions, and differential aperture photometry. 

A first full transit was observed on 2020 February 22 in the i’ 
filter with a sequence of 30s exposure time. As seen in the 
lower panel of Figure 1, TOI-1452 and TIC 420112587 were 
partially resolved. Using a circular aperture of 7”9 containing 
both stars, the transit was detected 53 min earlier than predicted 
by the TTF (2.20 early), causing us to miss observing a proper 
pre-ingress baseline (see Figure 3, upper left). Additional TESS 
data later confirmed that the period was slightly overestimated 
by SPOC (sectors 14-16 only) at the time of observations, 
explaining why the transit arrived “early”. The transit timing 
was also later confirmed by TESS sector 22 data, which was 
contemporaneous to this data set. Even without a proper 
baseline, this transit was particularly valuable because it 
allowed us to reject the NEB false-positive scenario and to 
determine, using point-spread function (PSF) fitting (see 
Section 4.1), that the signal originated from TOI-1452. 

A second full transit of TOI-1452.01 was observed on 2021 
March 4 in i’, using a 10s exposure time sequence. With a 
combined aperture of 8/4, the transit was detected in time 
according to the TTF. In addition to a standard airmass linear 
detrending (also performed for the first transit), we used the 
Width (mean of the x- and y-direction FWHM) detrending 
option in AI J. This was necessary to account for flux loss when 
the seeing was worse for certain exposures in the sequence, 
without increasing the aperture radius and dealing with flux 
contamination from a third star. This made sure that the transit 
depth was consistent with the one derived from the first PESTO 
observation, when the overall seeing was better (see Table 2). 

The OMM-PESTO observations are summarized in Table 2. 
The resulting aperture photometry transits are shown in 
Figure 3 and were included with the 32 TESS transits in our 
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Figure 3. Ground-based transit monitoring of TOI-1452.01 with the PESTO camera installed at OMM. Left panels: Individual transit observations in the i’ band on 
2020 February 22 UT and 2021 March 4 UT. Right panel: Transits from the left panels folded in phase, with black points representing the binned photometry 
(8 minute phase bin). The red curve shows the best-fit transit model (described in Section 4.3), with the 68% confidence interval envelope in transparent red. The 


residuals are shown below the phase-folded transit. 


joint analysis (transit and RV data sets) presented in 
Section 4.3. 


2.3. MuSCATS3 Transit Monitoring 


A full transit of TOI-1452.01 was observed on 2021 
September 8 with the multiband imager MuSCAT3 (Narita 
et al. 2020) on the 2m Faulkes Telescope North (FTN) of Las 
Cumbres Observatory (LCO) at Haleakala observatory, Hawaii. 
MuSCATS3 has four optical channels, each of which is equipped 
with a 2 kx2 k CCD camera with a pixel scale of 0266 pixel ', 
enabling g’-, r’-, i’-, and z,-band simultaneous imaging. For 
transit monitoring, significant chromaticity in the transit depths 
could indicate a diluted eclipsing binary. The exposure times 
were set at 35, 12, 6, and 3s for the g’, r’, i’, and z, bands, 
respectively. The observations were performed in focus to 
spatially resolve the host star from the nearby companion star at 
32, resulting in the FWHM of stellar point-spread function of 
3-8 pixels (0” 8-2”0) depending on the airmass and band (see 
Table 2). 

The obtained images were calibrated by the BANZAI 
pipeline (McCully et al. 2018). We performed aperture 
photometry on the calibrated images using a custom pipeline 
(Fukui et al. 2011) with aperture radii of 5 pixels, or 1/3, for all 
bands, which is almost free from flux contamination from the 
nearby companion star. For each band, we extracted the light 
curve using different sets of comparison stars but have found 
that using only the companion (TIC 420112587) as a reference 
produced the minimum point-to-point dispersion. As both the 
companion and the target stars have a similar color and are 
close to each other, the attenuation by the atmosphere is almost 
identical, so that we can safely assume that any second-order 
extinction effect (airmass dependent) is almost negligible. We 
clearly detected the transit on the target star in all bands, as 
shown in Figure 4, providing further unambiguous evidence 


Table 2 
Summary of the Ground-based Transit Monitoring of TOI-1452 

UT Date Camera Filter PSF Aperture 

FWHM (”) Size () 

OMM 1.6m 
2020-02-22 PESTO i 2.9 7.9" 
2021-03-04 PESTO i! 3.9 8.4" 
LCO-FTN 2m 

2021-09-08 MuSCAT3 g! 1.8 1.3 
2021-09-08 MuSCAT3 r! 1,3 1.3 
2021-09-08 MuSCAT3 i! 1.3 1.3 
2021-09-08 MuSCAT3 2s 1.0 1.3 
Note. 


* Using an aperture containing TOI-1452 and TIC 420112587. 


that TOI-1452 hosts a transiting object. A summary of this data 
set is provided in Table 2. 


2.4. Keck II/NIRC2 High-resolution Imaging 


One or more unresolved sources not in Gaia EDR3 could 
still be located close to TOI-1452, whether gravitationally 
bound or not. A blended eclipsing binary (BEB) could indicate 
a false-positive detection, and any other flux source would lead 
to underestimate the size of the transiting object in the TESS, 
PESTO, and MuSCAT3 light curves. For these reasons, we 
searched for subarcsecond sources around TOI-1452 with the 
NIRC2 adaptive optics imaging camera installed on the 10m 
Keck II telescope. The images were acquired on 2020 May 28 
in the K band with a spatial resolution of 0701 per pixel, 
integration time per coadd of 1.6s, mean PSF FWHM of 
0”061, and airmass of 1.69. Figure 5 shows the 5c contrast 
curve of TOI-1452, revealing that no additional companion is 
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Figure 4. Ground-based transit follow-up of TOI-1452.01 on 2021 September 8 UT with the multifilter (g’, r’, i’, z.) MuSCAT3 instrument installed on LCO-FTN at 
Haleakala Observatory. For each corresponding filter, the black points depict the binned photometry (8 minute temporal bin). The color coded curves correspond to 
each filter’s best-fit transit model (described in Section 4.3), with their respective 68% confidence interval envelope in lighter shade. The residuals are shown below 


each phase-folded transit. 


detected with a contrast ratio AK < 5.429 for separation greater 
than 0”5. Although in the K band, this contrast limit is similar 
to the difference in magnitude required (AT = 5.55) for a 50% 
depth BEB to mimic a 3 ppt transit in the TESS light curve. 
Following the procedure of Lillo-Box et al. (2014), we 
calculated the probability of contamination from a blended 
source due to a random alignment inside 0”5. For this, we 
simulated the galactic stellar population in a region near the 
target with TRILEGAL (Girardi et al. 2012), using their default 
bulge, halo, disk (thin and thick) parameters and the log-normal 
initial mass function of Chabrier (2001). The probability of an 
undetected source with AK <5.55 inside 0”5 is less than 
0.04%, so we can safely assume that the transit signal is not 
produced by a BEB or significantly diluted by a _back- 
ground star. 


2.5. SPIRou Velocimetry 


TOI-1452 was observed at 53 epochs from 2020 June 4 to 
2020 October 8 with the near-infrared (0.98-2.5 wm) SPIRou 
spectropolarimeter (Donati et al. 2018, 2020) mounted on the 
3.6 m Canada—France—Hawaii telescope (CFHT). The observa- 
tions were conducted as part of the ongoing SPIRou Legacy 
Survey (SLS; Donati et al. 2020), more precisely its Transit 
Follow-up program (SLS-WP2), which aims to characterize 
exoplanets orbiting low-mass stars revealed by photometric 
surveys such as TESS. SLS-WP2 has thus far allowed the 
characterization of the brown dwarf TOI-1278 B (Artigau et al. 
2021), the sub-Neptune TOIJ-1759b (Martioli et al. 2022) 
transiting M dwarfs, and the studies of the transiting planets 
HD 189733b (Moutou et al. 2020; Boucher et al. 2021) and AU 
Mic b (Martioli et al. 2020). 

SPIRou offers simultaneous high-resolution spectroscopy 
and polarimetry, with a spectral resolving power R ~ 70,000. 
Each epoch measurement consisted of four consecutive 15 
minutes exposures, i.e., a polarimetric sequence, with two 
rotating Fresnel rhombs varying positions between the 
exposures. During such a sequence, the two science fibers, A 
and B, each receive orthogonal polarization states, giving 
access to the circular polarization and total intensity of the light 
beam (Stokes V and J). A total of 212 spectra were collected, 
with SNR per spectral element (~2.2 km s | pixel! for 
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Figure 5. K-band 5o contrast curve of TOI-1452 from Keck II/NIRC2 
adaptive optics imaging. No close companion is detected. 


SPIRou) between 20 and 65 (median of 55) near 1.6 zm. Four 
individual spectra were rejected; one due to loss of guiding, 
three others because of high extinction (clouds). A single 
polarimetric sequence of four 15 minute exposures was also 
acquired on the 3”2 companion (TIC 420112587) on 2021 
April 22, principally to check its rotation profile and magnetic 
activity level (see Section 3.5). The diameter of the SPIRou 
fiber is 1!’33 and the typical seeing during the observations of 
TOI-1452 was 0”8. We measure no correlation between the 
radial velocity residuals (Keplerian and activity models 
described in Section 4.3) and the seeing, suggesting that any 
effect of contamination from the companion was negligible. 
Per-epoch RV measurements consisted of taking the error 
weighted mean of the individual observations within a 
polarimetric sequence. The data analysis presented in 
Section 4.3 was performed on the unbinned RVs, but we show 
the per-epoch average to facilitate visualization. 

The SPIRou data were reduced with APERO v0.7.194 
(N. J. Cook et al. 2022, in preparation). In brief, APERO starts 
by correcting known H4RG infrared detector defects (Artigau 
et al. 2018), then proceeds to identify bad pixels, locate each 
spectral order on the image, calculate the shape of the 
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instrument pupil slicer (Micheau et al. 2018), and finally 
determine from nightly calibration sequences the flat and blaze 
corrections to apply. Once this preprocessing and calibration 
step is completed, APERO performs an optimal flux extraction 
(Horne 1986) in both science channels, separately (fibers A and 
B) and together (AB), as well as in the simultaneous calibration 
channel (fiber C). The extracted 2D images (49 orders by 4088 
pixels) are then spectral flat fielded, as well as corrected for 
thermal background and for any leakage from the calibration 
channel to the science ones. A nightly pixel-to-wavelength 
solution is applied using a combination of a UNe hollow- 
cathode lamp and a Fabry—Perot (FP), as described in Hobson 
et al. (2021). APERO uses the simultaneous FP measurements 
from fiber C to calculate drifts between individual science 
frames relative to the nightly wavelength solution (typically 
below 2m s_'). Finally, a telluric absorption and night-sky 
emission correction is applied in a two-step process. The 
science frames are first pre-cleaned with a TAPAS (Bertaux 
et al. 2014) absorption model that leaves percent-level residuals 
for deep (>50%) H2O and dry absorption features (e.g., CHy, 
Ox, CO2, NO, and O3). Then, a telluric residuals model with 3 
degrees of freedom per pixel (optical depths for the H,O and 
dry components and a constant) is fitted to the pre-cleaned 
spectra. The grid of telluric models was generated from a set of 
rapidly rotating hot stars observed with SPIRou at various 
airmass, water columns, and dry absorptions, producing 
telluric-corrected spectra with final residuals at the level of 
the PCA-based method of Artigau et al. (2014). 

Radial velocity measurements were obtained from the telluric- 
corrected spectra using the novel line-by-line (LBL) method 
(Artigau et al. 2022). The LBL formalism is based on the 
Bouchy et al. (2001) framework, in which Doppler shifts are 
inferred for individual spectral lines (~16,000 for an M dwarf 
observed with SPIRou) as opposed to a given spectral range. As 
in Bouchy et al. (2001), such calculations require a noiseless 
template as velocities are derived from the comparison between 
the residuals (observed spectrum minus template) and the 
derivative of the template. For a given observed star, one uses in 
practice a high SNR combined spectrum as a template, so that 
any remaining noise is small compared to that of an individual 
spectrum. For TOI-1452, the combined spectrum produced by 
APERO did not reach a SNR as high as other bright standard 
stars observed in the SLS. Moreover, TOI-1452 is located near 
the North ecliptic pole, meaning its yearly barycentric Earth 
radial velocity (BERV) variation is small. Our observations with 
SPIRou covered BERV excursions between 1.7 and 4.8 km a, 
which is not ideal to filter out telluric lines (i.e., stellar lines do 
not move a lot with respect to the telluric lines), producing a 
template that still contains some telluric artifacts. For these 
reasons, we used the template of Gl 699, a standard star 
monitored with SPIRou for 2.5 yr with a spectral type (M4V) 
similar to that of TOI-1452 (M4 + 0.5, see Section 3.1) and a 
good BERV coverage (426kms°'). 

For each spectrum, the LBL algorithm combines thousands 
of per-line velocities into a single RV measurement, with per- 
line uncertainties varying from 50m s_' for the strongest 
features to tens of kms_' for the shallow ones. This is achieved 
using a simple mixture model: per-line velocities either 
originate from a Gaussian distribution centered on the mean 
velocity, with a standard deviation derived from Bouchy et al. 
(2001), or they arise from another distribution, namely, that of 
high-sigma outliers, whose plausible causes are diverse 
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(persisting bad pixels, cosmic rays, telluric residuals, etc.). 
Lastly, the LBL RVs are corrected for the instrumental day-to- 
day drift measured by the FP and for a long-term zero-point 
obtained with a Gaussian process regression using the most 
observed stars in the SLS. This zero-point calibration is similar 
to Courcol et al. (2015) for the SOPHIE spectrograph but will be 
described in more details in a forthcoming publication 
(T. Vandal et al. 2022, in preparation). The comparison between 
the LBL and other methods such as the cross-correlation 
function and template matching is discussed in Martioli et al. 
(2022) and in Artigau et al. (2022). The final SPIRou radial 
velocities of TOIJ-1452 are listed in Table D1, with typical 
precision of 8.0ms_' per exposure, or 4.0ms' per epoch. 


2.6. IRD Velocimetry 


Seven high-resolution spectra of TOI-1452 were obtained 
with the InfraRed Doppler (IRD) spectrograph on the Subaru 
8.2m telescope (Tamura et al. 2012; Kotani et al. 2018) 
between 2020 September 26 and 2021 June 25. IRD covers the 
near-infrared wavelengths between 970 nm and 1730 nm, with 
a spectral resolution R ~ 70, 000. For accurate RV measure- 
ments, stellar spectra were obtained simultaneously with the 
reference spectra of the laser-frequency comb (LFC). The 
integration times were set to 600-1500 sec, depending on the 
available observing time slots and sky conditions. The IRD 
fiber has a 0/48 diameter, so that flux contamination from the 
companion star is not an issue. 

The raw IRD data were reduced following the standard 
procedure of Hirano et al. (2020). We extracted wavelength- 
calibrated one-dimensional spectra for TOI-1452, as well as for 
the simultaneously injected LFC. The typical SNR of the TOI- 
1452 extracted spectra was 60-70 per pixel around 1000 nm. 
To measure precise RVs for TOI-1452, the reduced spectra 
were put into the RV analysis pipeline for IRD (Hirano et al. 
2020). This pipeline fits each small spectral segment of the 
observed spectra by the forward-modeling technique, taking 
into account the instantaneous variations of Earth’s atmo- 
spheric features as well as the instrumental profile of the 
spectrograph (which is estimated based on each laser-comb 
spectrum). The seven IRD RV measurements have an overall 
precision of 4.0m s’' and are given in Table D1. 


3. Stellar Characterization 
3.1. TOI-1452 (TIC 420112589) 


The star TOI-1452 (TIC 420112589) is a nearby M dwarf at 
a distance of 30.504+ 0.013 pe (Gaia Collaboration et al. 
2021). This star does not belong to any known young stellar 
moving groups, with a very high probability (>99.9%) of being 
a field star (Gagné et al. 2018). The presence of flares and 
short-period sinusoidal signal in the TESS PDCSAP data (see 
Figures 2 and A1) cannot be attributed with certainty to TOI- 
1452, due to flux contamination from multiple nearby objects. 
An analysis of the photometric variations is presented in 
Section 4.2, but we note that the polarimetric data from SPIRou 
reveal no important surface magnetic field variations (see 
Section 3.5), suggesting that TOI-1452 is relatively quiet, with 
a rotation period probably much longer than the modulation 
seen in the TESS light curve. 

As discussed in Section 3.4, we measure an effective 
temperature of 3185 +50K for TOI-1452 using the SPIRou 
combined spectrum, from which a spectral type (SpT) between 
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M4 and M4-.5 is inferred based on Table 5 of Pecaut & Mamajek 
(2013). We also considered the Gaia DR2 color to SpT relation of 
Kiman et al. (2019), more specifically the G— Gpp relationship, for 
an independent SpT determination. From this relationship, the Gaia 
magnitudes and their respective uncertainties, we obtain a SpT of 
M3.7 + 0.6. The same M4 spectral type was derived through a 
visual comparison of the SPIRou combined spectrum, degraded to 
a lower resolution (R ~ 5 000), with spectral type standards of the 
IRTF spectra library (Cushing et al. 2005; Rayner et al. 2009). 
Considering all these estimates, we adopt a SpT of M4 + 0.5. 

The mass of TOI-1452 was inferred from the Mann et al. 
(2019) absolute K, magnitude (Mx,) to M, relation for M 
dwarfs. Taking into consideration the dispersion of this 
relation, the K, magnitude, the distance, and their corresp- 
onding uncertainties, a mass of 0.249 + 0.008 Mq is obtained. 
A similar approach was used to measure the stellar radius, this 
time using the Mx —R, relationship of Mann et al. (2015), from 
which we derive R,=0.275+0.009 Rs. Other physical 
parameters such as the surface gravity (log g), the mean density 
(p,), and the luminosity (L,) were determined from the M,, R,, 
and Tere estimates. The stellar parameters of TOI-1452 are 
summarized in Table 3. 


3.2. Bound Companion (TIC 420112587) 


TOI-1452 has a resolved companion (TIC 420112587) with 
several comparable photometric and astrometric measurements. 
The two objects have similar Gaia EDR3 magnitudes of 
G= 13.598 + 0.003 and G= 13.830+0.003 for TOI-1452 
and TIC 420112587, respectively. Their Gaia EDR3 parallaxes 
are identical (within the errors), 32.782 + 0.014 mas for TOI- 
1452 and 32.791 + 0.014 mas for TIC 420112587, indicating a 
very similar distance to these stars. Their projected angular 
separation is 3”182, which corresponds to a projected physical 
separation of ~97 au, using a common approximate distance of 
30.5 pe. The proper motion of TIC 420112587 is similar to that 
of TOI-1452, with ju, cos 6 and jzs5 within 15% for the two stars 
(see Tables 3 and 4). TOI-1452 and TIC 420112587 most 
likely form a visual binary, i.e., a resolved gravitationally 
bound system, which was previously reported in the TOI 
visual-binary catalog of Mugrauer & Michel (2020), as well as 
in the binary catalog based on Gaia EDR3 of El-Badry et al. 
(2021). From our single visit on the companion star with 
SPIRou, we measure an RV offset between TOI-1452 and TIC 
420112587 of —6.9kms~'. The SPIRou template spectrum of 
TIC 420112587 combining only four individual spectra at the 
same epoch and BERV does not allow for a similar spectral 
analysis as the one presented for TOI-1452 in Section 3.4. 

Using the empirical relationships of Mann et al. 
(2015, 2019), we obtain a mass of 0.226+ 0.006 M, and a 
radius of 0.254 + 0.008 R. for TIC 420112587. The mass ratio 
of the binary system is close to unity (¢ = 0.91 + 0.04), with 
TOI-1452 as the primary member. The projected physical 
separation and masses of the two stars imply an orbital period 
of about 1400 yr. The radial velocity variation expected from 
such an orbital motion and for a circular orbit is under 
~1.5ms_' over the span of our SPIRou RV observations. 
According to Table 5 of Pecaut & Mamajek (2013) and a 
T¢ = 3060 + 50 K derived from the spectral energy distribu- 
tion (see analysis below), TIC 420112587 has an M5 spectral 
type. The Kiman et al. (2019) G— Ggp relationship yields an 
SpT of M4.0 + 0.6, so we adopt an intermediate spectral type 
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Table 3 
TOI-1452 Stellar Properties 

Parameter Value Ref. 

Designations 
TIC 420112589 1 
TOI 1452 1 
2MASS J19204172+7311434 2 
UCAC4 816-023943 3 
Gaia EDR3 2264839957 167921024 4 

Astrometry 
RA (J2016.0) 19:20:41.75 4 
DEC (J2016.0) +73:11:42.35 4 
H, cos 6 (mas/yr) 7.800 + 0.017 4 
fis (mas/yr) —74.076 + 0.017 4 
™ (mas) 32.7823 + 0.0140 4 
Distance (pc) 30.5043 + 0.0130 4 

Stellar parameters 

Tere (K) 3185 + 50 5 
SpT M4 + 0.5 5 
[M/H] —0.07 + 0.02 =) 
M,. (Mo) 0.249 + 0.008 a 
R, (Ro) 0.275 + 0.009 5 
log g (dex) 4.95 + 0.03 5 
p, (g em~4) 16.8+1.9 5 
L, (Ls) 0.0070 + 0.0006 5 

Photometry 
B 15.94 + 0.03 1 
Vv 14.35 + 0.12 1 
Gpp 15.222 + 0.004 4 
G 13.598 + 0.003 4 
Grp 12.362 + 0.004 4 
T 12.295 + 0.007 1 
g 15.580 + 0.002 6 
r 14.383 + 0.007 6 
i 12.873" 6 
Zz 12.272°* 6 
y 11.875 + 0.020 6 
J 10.604 + 0.058 2 
H 10.026 + 0.058 2 
K, 9.740 + 0.046 2 
wl 8.938 + 0.023” 7 
Ww2 8.760 + 0.019° 7 
W3 8.686 + 0.023” 7 
w4 8.46 + 0.29° 7 
Notes. 


*“ The uncertainty was not indicated. 

> WISE magnitudes include the flux from TOI-1452 and TIC 420112587. 
References: (1) TIC (Stassun et al. 2019). (2) Two Micron All Sky Survey 
(2MASS; Skrutskie et al. 2006). (3) UCAC4 (Zacharias et al. 2013). (4) Gaia 
EDR3 (Gaia Collaboration et al. 2021). (5) This work. (6) Pan-STARRS1 DR2 
(Chambers et al. 2016). (7) AIIWISE (Wright et al. 2010). 


of M4.5+0.5. A summary of the stellar properties of TIC 
420112587 is presented in Table 4. 


3.3. Spectral Energy Distribution Fit 


As an independent determination of the basic stellar 
parameters, as well as to estimate the contaminating flux from 
the nearby companion star, we performed an analysis of the 
broadband spectral energy distribution (SED) of the stars 
together with the Gaia EDR3 parallax (with no systematic 
offset applied; e.g., Stassun & Torres 2021), following the 
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Table 4 
TIC 420112587 Stellar Properties 

Parameter Value Ref. 

Designations 
TIC 420112587 1 
TOI 1760 1 
2MASS J19204172+7311467 2 
Gaia EDR3 2264839952875245696 3 

Astrometry 
R.A. (J2016.0) 19:20:41.76 3 
Decl. (J2016.0) +73:11:45.53 3 
Hi, cos 6 (mas/yr) 6.845 + 0.017 3 
Js (mas/yr) —82.216 + 0.017 3 
m™ (mas) 32.7913 + 0.0141 3 
Distance (pc) 30.4959 + 0.0131 3 

Stellar parameters 

Tepe (RK) 3060 + 50 4 
SpT M4.5 + 0.5 4 
M, (M5) 0.226 + 0.006 4 
R, (Ro) 0.254 + 0.008 4 
log g (dex) 4.98 + 0.03 + 
ps (g cm” *) 195+18 4 
L, (Le) 0.0051 + 0.0005 4 

Photometry 
B 15.76 + 0.17 1 
Vv 13.99 + 0.2 1 
Gpp 15.512 + 0.005 3 
G 13.830 + 0.003 3 
Grp 12.576 + 0.004 3 
T 12.499 + 0.008 1 
g 15.890 + 0.002 5 
r 14.659 + 0.003 5 
i 13.153 + 0.002 5 
Zz 12.456 + 0.021 5 
y 12.111 + 0.007 5 
J 10.795 + 0.027 2 
H 10.257 + 0.031 2 
Ks 9.944 + 0.023 2 
Note. 


References: (1) TIC (Stassun et al. 2019). (2) 2MASS (Skrutskie et al. 2006). 
(3) Gaia EDR3 (Gaia Collaboration et al. 2021). (4) This work. (5) Pan- 
STARRS1 DR2 (Chambers et al. 2016). 


procedures described in Stassun & Torres (2016) and Stassun 
et al. (2017, 2018a). For both stars, we pulled the JHKs 
magnitudes from 2MASS, the W1—W4 magnitudes from Wide- 
field Infrared Survey Explorer (WISE), and the grizy 
magnitudes from the Panoramic Survey Telescope and Rapid 
Response System (Pan-STARRS). Together, the available 
photometry spans the full stellar SED over the wavelength 
range 0.4-10 yum (see Figure 6). We excluded the WISE 
photometry from the initial fitting because the two stars are 
blended in WISE, such that the catalog photometry in fact 
represents the sum of the fluxes of both stars. 

For each star, we performed a fit using NextGen stellar 
atmosphere models (Hauschildt et al. 1999), with the free 
parameters being the effective temperature (T.,¢) and metallicity 
({Fe/H]). The remaining free parameter is the extinction Ay, 
which we fixed at zero due to the stars’ proximity. The 
resulting fit for TOI-1452 (Figure 6) has a reduced 7 of 1.8 
with Ter = 3100 +50 K and [Fe/H] =0.0+0.5. Integrating 
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Figure 6. SEDs of TOI-1452 and TIC 420112587 fitted with a NextGen 
atmospheric model, respectively depicted with a black and a cyan line. 
Photometric measurements are represented with red (TOI-1452) and blue (TIC 
420112587) error bars, where the horizontal bars represent the effective widths 
of the bandpasses and the small blue dots represent the model fluxes of TOI- 
1452 for comparison to the red symbols. The WISE measurements (3.4, 4.6, 
12, and 22 yum), excluded from this analysis for including both star fluxes, 
agree with the sum of the two SEDs (green). 


the model SED gives the bolometric flux at Earth, 
Fyo1 = 2.34 £0.11 x 1071! ergs-'cm~?. Taking the F,,) and 
Tee together with the Gaia parallax, gives the stellar radius, 
R,, = 0.286 + 0.011 R,. This independent radius measurement 
is consistent, although slightly less precise, with the one 
derived using Mann et al. (2015; R, = 0.275 + 0.009 R35). 

Similarly, the resulting parameters for the companion star 
from the SED fit are Tos = 3060 + 50 K, [Fe/H] =0.0+0.5, 
and R, = 0.263 + 0.010 R.. This radius estimate is again fully 
consistent with the value derived from empirical relation 
(R,, = 0.254 + 0.008 Ro). The sum of the two stellar models is 
compared to the combined WISE fluxes in Figure 6, showing 
good agreement. Integrating the companion SED within the 
TESS bandpass yields a flux ratio (companion relative to TOI- 
1452) of 0.77 + 0.03. Note that the flux ratio derived strictly 
from the T magnitudes from the TIC is 0.829 + 0.002. In the 
event that the PDCSAP overestimated the dilution correction 
for TIC 420112587, this difference in flux ratio would imply a 
~1.7% overestimation of the planetary radius. 


3.4. Stellar Parameters from SPIRou Spectra 


The high-resolution combined spectrum of TOI-1452 from 
SPIRou lets us determine T.¢¢ and the abundance of several 
elements with relatively good accuracy. This work follows the 
methodology of F. Jahandar et al. (2022, in preparation), 
which we briefly summarize here. Because models and 
observations can often show significant discrepancies (e.g., 
continuum mismatch in the Y and J bands), we only select for 
the fitting analysis a subset of relatively strong lines that are 
matching the models. The selected lines are then divided into 
several groups of 15 lines, each analyzed independently 
through a chi-squared fitting routine to infer both Trp and [M/ 
H] for all groups. The spectrum is compared with a grid of 
ACES models (Allard et al. 2012a; Husser et al. 2013). The 
advantage of this method is that it yields several (typically 15) 
independent measurements that can be used to characterize 
the inherent uncertainties associated with the fitting 
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Figure 7. SPIRou observation of the Al I line (1675.514 nm) of TOI-1452 
(black dashed line). Top panel: The solid lines represent the ACES models for a 
fixed metallicity of 0.1 dex and T.4 values of 3060 K, 3160 K, and 3260 K. 
Bottom panel: Same as top panel, but the ACES models have a fixed 
Tere = 3160 K and metallicity values of —0.2 dex, 0.1 dex and 0.5 dex. These 
plots illustrate the good sensitivity of near-IR high-resolution spectroscopy for 
constraining both the metallicity and effective temperature of M dwarfs. 
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procedure. This analysis applied to the TOI-1452 spectrum 
yields Terp= 3185450 K and [M/H] = —0.07 + 0.02, in 
good agreement with the parameters derived from the SED 
fitting analysis. The quoted uncertainty for Ter is internal to 
our fitting methodology and ignore potential systematic 
differences with bolometric T.- estimates based on interfero- 
metric measurements. While our Trp estimates have yet to be 
calibrated with bolometric Tyr, it is empirically demonstrated 
that 50-60 K is a typical uncertainty derived from atmosphere 
models (e.g., Mann et al. 2013, 2015). We thus adopt 50 K for 
our Tz uncertainty, a conservative value given that the 
temperature derived from the SPIRou spectrum is inferred 
from several independent measurements. An illustration of 
the temperature and abundance sensitivity for an Al I line 
(at 1675.514nm) is shown in Figure 7. In practice, several 
tens of lines are used to derive Terr. 

Once Tere is determined, one can then proceed, through a 
similar procedure, to determine the abundance of all individual 
lines of a given element. The high-resolution near-infrared 
spectrum of an M dwarf is characterized by several hundreds of 
relatively strong OH lines. By selecting only those that are well 
isolated, i.e., with no known spectral features within a few 
pixels using the PHOENIX/BT-Settl (Allard et al. 
2012b, 2013) and the National Institute of Standards and 
Technology (Kramida et al. 2021) line lists, we find 72 OH 
lines, whose individual abundance can be used to quantify the 
inherent, per-line uncertainty of this method. This uncertainty 
obviously does not consider any possible systematic errors 
associated with the ACES atmosphere models. The 72 
independent OH abundance measurements are presented in 
Figure 8, showing a good match with a Gaussian distribution 
with standard deviation 0.13dex. For all elements and 
molecules detected in TOI-1452, we list the average abundance 
of all lines in Table 5 (see also Figure 9). For chemical species 
with only one line, we adopt an uncertainty of 0.13 dex from 
the OH distribution. We report abundances for Fe, Mg, and Si 
that constitute the bulk material of an exoplanet core and 
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Figure 8. Distribution of OH abundance measurements from 72 isolated OH 
lines in the SPIRou combined spectrum of TOI-1452. The inherent, per-line 
uncertainty (0.13 dex) is inferred from a Gaussian fit of this distribution. This 
uncertainty does not take into account possible systematic errors of the stellar 
atmosphere models. 


Table 5 
Stellar Abundance of TOI-1452 for Various Chemical Species Measured by 
SPIRou 
Element [X/H] o # of Lines 
Fe I —0.07 0.03 38 
All 0.07 0.11 4 
Mg I 0.02 0.07 a 
Sil 0.11 0.13 1 
Til —0.31 0.06 10 
Ca I 0.01 0.12 2 
Cri 0.04 0.06 4 
KI 0.03 0.13 1 
or —0.24 0.02 72 
CI —0.17 0.05 11 
NI —0.12 0.13 1 
NaI —0.22 0.06 
Z>" —0.07 0.02 
Notes. 


* The oxygen abundance is inferred from OH lines. 
2 Average abundance of all elements. 


mantle. The overall metallicity ((M/H]) and its corresponding 
error are determined by averaging the final abundance of each 
element in Table 5, assuming a common uncertainty for all 
elements taken as the median of all individual uncertainties. 
This approach is chosen to avoid putting too much weight on 
the oxygen abundance characterized by a small uncertainty. 


3.5. Spectropolarimetry with SPIRou 


The combination of the four exposures within a polarimetric 
sequence obtained with SPIRou yields the circular polarization 
profile at the surface of the star (Donati et al. 2020). The 
intensity (Stokes J), circular (Stokes V), and null polarization 
spectra were generated in APERO following the spirou- 
polarimetry code.** We applied the least-square deconvo- 
lution (LSD) method of Donati et al. (1997), also outlined in 
Martioli et al. (2020), to compute the average I and V profiles. 


ae github.com/edermartioli/spirou-polarimetry 
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Figure 9. Chemical abundances of 12 different elements via line-by-line 
spectroscopy of 151 spectral features. The dashed line represents the average 
metallicity of the star corresponding to the average of all elements. 


We used the Vienna Atomic Line Database (Piskunov et al. 
1995) and a MARCS atmosphere model (Gustafsson et al. 2008) 
with Te = 3000 K and log g = 5.0 dex to search for valid atomic 
features. Lines deeper than 3% and with a known Landé factor 
were selected to produce a line mask of 955 atomic lines, used in 
this LSD analysis of TOI-1452. An estimate of the longitudinal 
magnetic field (B,) at the stellar surface can then be obtained 
using Equation (5) of Donati et al. (1997), combining the Stokes 
I and V LSD profiles, the mean Landé factor of 1.24, and the 
mean wavelength of 1604.59 nm. By doing this over multiple 
epochs, one can monitor the large-scale surface magnetic field, 
expected to vary with the rotation of the star. The polarimetric 
capabilities of SPIRou can thus serve as a useful activity tracer 
simultaneous to the RV measurements, as demonstrated in 
Martioli et al. (2022), where the rotation period of the 
moderately active MO star TOI-1759 (Rot = 35.65" 912 days) 
was determined from the By, time series. We obtained 
independent and consistent values for the By of TOI-1452 using 
the Libre—Esprit pipeline (Donati et al. 1997, 2020) but present 
below the values from the APERO pipeline. 

The B;, time series of TOI-1452 is presented in Figure 10. A 
simple Lomb-Scargle periodogram analysis shows no obvious 
periodicity. The By data do not favor a sinusoidal model, which 
could be associated with stellar rotation, over a constant 
magnetic field (mean By =—3.8 + 1.8 G). The small variation 
of By suggests that the field is intrinsically weak (quiet star), or 
that it is strongly axisymmetric with respect to the rotation axis. 
Alternatively, the rotation period of TOI-1452 could be longer 
than the 4 month span of our observations, but this is close to 
the largest known period for M dwarfs (Newton et al. 2018). 

This LSD analysis was also applied on the companion star TIC 
420112587 using the single polarimetric sequence acquired with 
SPIRou. We report a polarimetric signal and a By consistent with a 
null value, indicating that the companion is also probably inactive. 


4. Data Analysis & Results 
4.1. Determining the Transit Origin with PSF Photometry 


The objective of the first OMM-PESTO transit follow-up 
was to establish the origin of the TESS signal, particularly 
between the target (TOI-1452) and its companion (TIC 
420112587). Standard aperture photometry ruled out any 
NEB in the FOV but was unable to isolate the transit between 
the two stars, as they were only partially resolved. We therefore 
had to rely on a different method using PSF fitting to extract the 
relative flux of both stars. We used the photutils (Bradley 
et al. 2020) package to perform the DAOPHOT (Stetson 1987) 
PSF photometry algorithm. This was achieved by fitting the 
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PSFs with an effective PSF model (ePSF) generated in 
photutils using the six stars with the highest SNR in the 
FOV (excluding our targets), then integrating the best-fit 
models over pixels containing the stars’ signal. 

The intent here was not to produce a precise uncontaminated 
light curve but rather to detect any flux deficit (or excess) that 
would indicate from which star the transit originates. We thus 
inspected the TOI-1452 to TIC 420112587 flux ratio as a function 
of time, normalized to unity outside of transit. The resulting relative 
light curve is presented in Figure 11 and shows a flux deficit on 
TOI-1452 during transit. We did not fit a transit model on this light 
curve, as it is less precise than the one obtained using a combined 
circular aperture (Figure 3). We nonetheless measure a mean 
relative flux deficit of 2.33 + 0.43 ppt, which is an approximation 
of the uncontaminated transit depth. This flux deficit is comparable 
in amplitude to the diluted corrected TESS depth (3.31 + 0.19 ppt) 
and was detected with a confidence level sufficiently high (>5c) to 
conclude that TOI-1452 was the source of the transit and justify an 
RV monitoring campaign on this star, starting with SPIRou in June 
2020. Later, the MuSCAT3 photometry was able to resolve TOI- 
1452 and TIC 420112587 and unambiguously identify that the 
former star hosts a transiting object. 


4.2. TESS Light Curve Analysis 


The TOI-1452 PDCSAP light curve (Figures 2 and Al) 
features stellar flares with amplitude of a few percents and ppt- 
level sinusoidal variations. A strong peak at 0.93 days appears 
in the Lomb-Scargle periodogram of the multiyear light curve, 
as well as in all individual sectors. However, computing the 
autocorrelation function, which is more reliable for accurate 
photometric rotation period determination (McQuillan et al. 
2013), would often find a period of 1.9 days (2 x 0.93 days) 
depending on the sector. As the PDCSAP data are corrected for 
systematic trends, it is unlikely that such corrections sig- 
nificantly perturb those short-term flares and_ sinusoids. 
Regardless of the origin of these signals (TOI-1452, TIC 
420112587, or any contaminating star), it is crucial to remove 
the periodic variations to accurately measure the transit 
parameters. To accomplish this, we adopted a sequential 
approach where we first correct the PDCSAP data using a 
Gaussian process (GP), then fit the 32 corrected transits with a 
model. The details of the GP regression are presented below, 
while the transit modeling is described in Section 4.3. 

We started by masking the epochs of transit and removing 
outliers from the PDCSAP light curve with sigma clipping. It was 
determined that a 3.50 clipping was robust enough to remove both 
obvious outliers and stellar flares. This sigma clipping removed 
less than 0.2% of the out-of-transit data. Parts of sectors 21 and 47 
coinciding with TESS momentum dump events show large 
amplitude variations; those were considered to be nonastrophy- 
sical and were manually rejected. We also rejected data points in 
sectors 40 and 41, as they are isolated and have a median 
considerably different than unity. The data not considered in this 
analysis are either displayed in blue (transits) or in red (rejected) in 
Figures 2 and Al. The cleaned out-of-transit PDCSAP data set 
was too large (NV = 237 634) to be efficiently modeled with a GP. 
We therefore binned the data and instead used the corresponding 
1 hour effective cadence light curve (V = 7924). 

The GP regression was done with celerite2 (Foreman- 
Mackey et al. 2017; Foreman-Mackey 2018). We selected its 
RotationTerm kernel because it was specifically designed 
to model a range of quasiperiodic variability, from stellar 
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Figure 10. Top panel: Monitoring of the longitudinal magnetic field (By) of 
TOI-1452 with SPIRou. Bottom panel: Lomb-Scargle periodogram of the By 
time series. No clear periodic signal is detected, suggesting a relatively 
quiet star. 


rotation to pulsations. This kernel is the sum of two 
stochastically driven, damped harmonic oscillator (SHO) terms 
(SHOTerm) capturing both primary (Pcp) and secondary 
(Pgp/2) modes in Fourier space. The Fourier transform of 
the covariance function, known as the power spectral density 
(PSD), takes the following form: 


2 Si dae 
S(w) = j ee 
m (w? — wt)? + wi w?/Q? 


we Sw 
Nt w= w3)?? + w3 w?/Q?’ 


where each SHOTerm PSD is described by their respective 
power S,, Sy at w=0, their undamped angular frequency w), 
W , and their own quality factor Q);, Qo. 

As the periods of the two oscillators are separated by a factor 
of 2 (w2 = 2w,), the parameterization below reduces by one the 
number of free parameters: 


(1) 


a, = ¥S1w1Q1, (2) 

02 = JS2wW202, (3) 

ye (4) 
Wy 

nm = 282, (5) 
W2 

Pop = = =, 6) 

Wy WwW. 


where 0), 0 are the standard deviations (amplitudes) of the 
primary and secondary modes, 7;, 72 are the damping 
timescales of the primary and secondary oscillations, and Pgp 
is the undamped period of the primary mode. Note that these 
parameters differ slightly from the default RotationTerm 
kernel parameterization by making no assumptions on the 
relative amplitudes and quality factors between the two modes. 
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Figure 11. PSF photometry relative light curve (TOI-1452 to TIC 420112587 
flux ratio) from OMM-PESTO on 2020 February 22. The dotted black lines 
represent the ingress and egress of the transit estimated from contemporaneous 
TESS sector 22 data. The light curve is normalized with the out-of-transit 
median. A flux deficit comparable to the TESS transit depth is detected on TOI- 
1452 during transit. 


Our PDCSAP GP model consisted of the five hyperpara- 
meters above, plus an excess white noise term s. We sampled 
the posterior distributions of the parameters in their logarithmic 
form {In oj), Ingo, Int, Int, In Pop, Ins} using the Markov 
chain Monte Carlo (MCMC) package emcee (Foreman- 
Mackey et al. 2013) and a Bayesian formalism. We employed 
100 walkers and performed 100,000 steps with a burn-in of 
10,000. The number of steps was greater than 50 times the 
autocorrelation timescale for each parameter, which usually 
indicates a sufficient level of convergence (Sokal 1997; 
Foreman-Mackey et al. 2019). The adopted prior distributions 
and the posteriors median, 16th and 84th percentiles are 
reported in Table 6. The resulting mean GP prediction is shown 
in Figures 2 and Al superimposed on the original PDCSAP 
cadence. Even though the sinusoidal variations visually appear 
to repeat every ~0.93 day, our model converged to a very well 
constrained primary oscillation of 1.8680 + 0.0004 days, thus 
indicating significant power at the second harmonic. 

It is beyond the scope of this study to assess the exact cause 
of this strong and persistent signal, but we showed earlier that 
the SPIRou magnetic field constraints of TOI-1452 are 
inconsistent with a fast rotator and active object. Moreover, a 
1.9 days rotation period for TOI-1452 would correspond to a 
vsini of ~7kms_', readily detectable in the SPIRou 
combined spectrum. Instead, the mean line profile FWHM 
measured from the cross-correlation function (CCF) calculated 
in APERO suggests a slow rotator (ie., v sini < 2 kms_'). We 
repeated this step for the companion star from the single visit 
with SPIRou and also measured a FWHM consistent with 
vsini< 2kms_!. Thus, the rotation of the companion star is 
also most probably not causing this photometric signal. 


4.3. Joint Transit RV Fit 


In order to constrain the physical and orbital parameters of 
TOI-1452b, we conducted a joint analysis of the transits 
(TESS, OMM-PESTO, and MuSCAT3) and the RV data 
(SPIRou and IRD). The joint fit was performed using the 
juliet (Espinoza et al. 2019) package, which utilizes 
batman (Kreidberg 2015) to generate transit models and 
radvel (Fulton et al. 2018) to compute Keplerian RV models. 
The juliet framework implements nested sampling algo- 
rithms to sample posterior distributions, while also enabling 
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Table 6 
Prior and Posterior Distributions of the Quasiperiodic GP Model of the TOI- 
1452 TESS Light Curve (Details in Section 4.2) 


Parameter Prior Posterior 
Ing, U(—10, 0) —7.5103 
Ino U(—10, 0) —7.17$3 
In [7 /days] U(—10, 10) — 1.42733 
In [7 /days] U(—10, 10) TA 

In [Pop /days] U(—2, 5) 0.62487 0.000 
Ins U(—15, 0) =12,7248 


Note. U(a, b) is the uniform distribution between values a and b. 


model comparison via evaluations of the Bayesian log- 
evidence (InZ). We chose the dynesty (Speagle 2020) 
dynamic nested sampling option in juliet. Standard nested 
sampling (Skilling 2006) was designed to estimate evidences, 
not posteriors, and thus struggles with parameter estimation for 
complex distributions. Dynamic nested sampling (Higson et al. 
2019), on the other hand, adapts the number of live points 
based on the structure of the posteriors, providing parameter 
estimation comparable to MCMC algorithms. 

The transit and RV components of the joint fit have four 
parameters in common: the orbital period P, the time of inferior 
conjunction fo, the eccentricity e, and the argument of periastron 
w. For the transit modeling, we followed the parameterization 
from Espinoza (2018) of the impact parameter b and the planet- 
to-star radius ratio p= R,/R, to efficiently sample physically 
plausible values in the (b, p) space. Instead of fitting the scaled 
semimajor axis a/R,, we used the stellar density p, parameter- 
ization available in juliet. Fitting p, takes into account any 
prior information on the stellar mass and radius. We adopted a 
Gaussian prior on p, using the value and uncertainty in Table 3. 
Stellar limb-darkening effects in TESS, OMM-PESTO, and 
MuSCATS3 transits were modeled using the quadratic g, and q2 
parameters defined in Kipping (2013). For each instrument, we 
included in juliet a flux dilution factor D, a baseline flux M, 
and an extra jitter term 0. We set Drgss to 1 (no dilution), as the 
PDCSAP data are already corrected for crowding effects. The 
OMM.-PESTO light curve combines the flux of TOI-1452 and 
TIC 420112587, which requires an adequate Dpgsto factor to 
compensate for contamination. We thus constructed a Gaussian 
prior on Dpgsto With a mean value calculated with Equation (6) 
of Espinoza et al. (2019) and flux ratio derived from TOI-1452 
and TIC 420112587 magnitudes in the i band (from Tables 3 and 
4). The adopted prior on Dpgsto was (0.564, 0.0564?), that is, 
with a 10% standard deviation to account for errors on the 
magnitudes and deviations between i and i’. We also explored 
fixing Dpgsto to 0.564, while letting Drgss vary freely between 
0 and 2. Both approaches yielded a consistent measurement of 
the planetary radius (within 1c), indicating that the PDCSAP 
fluxes were in all likelihood properly corrected for contamina- 
tion. The dilution in the MuSCAT3 light curves was a priori 
unknown. However, it is expected that the g’ transit was more 
affected by dilution, as the seeing was worse for this filter (see 
Table 2). We adopted a conservative approach where a different 
Dywuscat3 is applied for each filter, with uniform priors between 
0.5 (twice the flux) and 1. 

The parameters specific to the RV Keplerian component 
were the semi-amplitude K, per-instrument offsets , and extra 
white noise terms o. We explored adding a global GP to model 
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common stellar activity signal in the SPIRou and IRD data. For 
this, we used the GP implementation in juliet that runs 
celerite (Foreman-Mackey et al. 2017). We chose the 
Matérn-3/2 approximation kernel, which takes the following 
form: 


kij(r) = AGI + 1/e)e2-9" (1 — 1fe)e@ +9], (7) 


where 7 = |t; — ¢)| is the time interval between data points i and 
J, Acp 1s the amplitude of the GP, w = J37/lgp, with lgp the 
timescale of the GP, and « is set to 0.01 (when «0, kj; 
converges to a Matém-3/2 kernel). We did not fit a per- 
instrument Acp and fcGp due to the limited number of RV 
measurements from IRD. We also considered choosing a 
quasiperiodic kernel in celerite instead (Equation (56) of 
Foreman-Mackey et al. 2017). As no clear periodicity was 
detected in the B, time series, or other activity indicators from 
the LBL such as the dLW metric (Zechmeister et al. 2018) or 
chromatic velocity slope changes, we applied a uniform prior 
on the stellar rotation period, namely, U/(0.1, 120) days. We 
found that the Matérn-3/2 kernel gave equivalent results with 
fewer hyperparameters needed (two instead of four) and that 
the quasiperiodic GP did not converge to a specific rotation 
period, showing no preference for a period of 0.93 days (or 2 
x0.93 days) as seen in TESS photometry. This is another 
indication that the sinusoidal signal in the out-of-transit 
PDCSAP data is probably not associated with TOI-1452 stellar 
activity. 

We examined the change in Bayesian log-evidence for a 
suite of joint models (M), all having an identical transit 
component. The “zero”-planet model (Mop) has a K fixed to 
Oms_', with only the RV offsets and extra white noise terms 
allowed to vary. This model tests whether the RV dispersion 
can be fully explained by white noise only, without questioning 
the transit detection. Single-planet models can either be with 
circular (Micp; e =0, w= 90°) or eccentric orbits (Miep; free 
e, w). Two additional models include a global RV activity GP 
(Mcp+cp and Miep+cp). To objectively assess the contribution 
from the IRD observations, we decided to apply this framework 
first on the SPIRou data individually and then using the full RV 
data set (SPIRou + IRD). 

For two competing models, the difference in log-evidence 
(A In Z) informs on the probability that one model supports the 
data better than the other. To interpret the significance of the 
A1nZ and select the “best” model, we followed the empirical 
scale introduced in Table 1 of Trotta (2008). A AlnZ > 5 
translates into “strong” evidence toward the model with the 
highest InZ. A 2.5 < AlnZ < 5 corresponds to “moderate” 
evidence, while AlnZ < 2.5 shows “weak” evidence at best, 
1.e., neither model should be favored in that case. 

Figure 12 shows the Bayesian log-evidence for different 
joint models and data sets. Note that the typical errors on the 
InZ computed by dynesty were 0.5, so that the AlnZ 
presented in Figure 12 have associated uncertainties of 0.7. We 
first observe that all planetary models are strongly favored 
(AlnZ > 5) compared to the “zero”-planet solution (Mop), 
providing quantitative evidence that the TOI-1452b Keplerian 
signal is detected in velocimetry, in phase with transit. There is 
also compelling evidence for models with an RV activity GP 
(Micp+cp and Mjep;cp), increasing InZ by approximately 10 
relative to Miep and Mjep. However, considering only the 
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Figure 12. Bayesian log-evidence (InZ) for different joint transit RV models 
(M) and RV data sets. The typical uncertainty on the A InZ is 0.7. The “zero”- 
planet model Mo, has a fixed K =0m sl Single-planet models Mj¢p and 
Mtep correspond to an RV component with a circular (e = 0, w= 90°) or 
eccentric (free e, w) orbit, respectively. Models Micpigp and Miep+cp add a 
Gaussian process to remove correlated noise in the RV data. 


SPIRou data yields similar or slightly larger AlnZ for all M 
compared to joint fits that include the seven IRD RVs. This 
suggests that the IRD observations do not significantly 
contribute to improve the Keplerian solution for TOI-1452b. 
The median RV uncertainty from IRD (4.03 ms_') is nearly 
identical to SPIRou (4.00 m s~'), but the point-to-point scatter 
(rms) is much larger: 12.71 ms! for IRD and 5.76ms_! for 
SPIRou. The models fail to capture the extra rms in the IRD 
data and are instead converging to solutions with white noise 
term comparable to the overall scatter (see Table B1). This is 
apparent in Figure 13 showing the RV component of the joint 
fit (model Mi.p+cp) using the full data set, with each 
instrument having their original error bar plotted. The IRD 
radial velocities were produced using the template spectrum of 
TOI-1452. As previously mentioned, this star has a small 
BERV excursion, which is not ideal for filtering out telluric 
lines. This may explain the increased dispersion in the resulting 
RVs, in this case, at a level much larger than the Keplerian 
signal. For this reason, we opted to present below the results 
using only the SPIRou RVs. We nonetheless provide all the 
relevant parameters of the RV modeling for the SPIRou only 
and SPIRou+IRD data sets in Table B1. 

The eccentric model Mjep,Gp produced the highest InZ 
(Figure 12) but with a Bayesian evidence indistinguishable 
from the circular model Mop+cp (AlnZ = 0.9). We report an 


eccentricity of 0.12*43, with e<0.32 at a 95% confidence, 
but argue that the simpler, circular model should be preferred at 
this point. The adopted priors and resulting posteriors of the 
Micp+cp fit are summarized in Table 7. The MuSCAT3 
photometric parameters are given in another table (Table 8) to 
facilitate comparison between filters. We measure dilution 
factors D for MuSCAT3 consistent with no dilution for the i’ 
transit, with a moderate level of contamination (~30%) in the 
g’ band. If we assume instead that the flux dilution was exactly 
zero for all filters, the uncorrected transit depths (uncorr.) 
presented in Table 8 show no sign of strong chromaticity. The 
best-fit transit models of the TESS, OMM-PESTO, and 
MuSCAT3 photometry are shown in Figures 2, 3, and 4, 
respectively. The best circular (Micp+cp) and eccentric 
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(Miep+cp) RV _ orbital fits of TOI-1452b are shown in 
Figure 14 in a phase-folded format. 


5. Discussion 
5.1. Planet Composition 


The analysis of the transit and RV data yields a mass of 
4.82+1.30M., and a radius of 1.672+0.071 Rg, which 
a convert into a planetary bulk density of 5. 6715 g 

a4 density similar to that of the Earth (5.5 gcm_ >) for a 
eike! that has more mass is indicative of an object composed 
of lighter material. We placed these measurements in a mass— 
radius diagram (Figure 15), along with various theoretical 
composition curves obtained by the interior structure model of 
Valencia et al. 2007, 2013; Plotnykov & Valencia 2020. We 
populated this diagram, and the figures in the discussion below, 
with data from the NASA Exoplanet Archive (Akeson et al. 
2013) using the exofile tool.*° 

By comparing the mass and size of TOI-1452b to theoretical 
mass-radius curves in Figure 15, we see three possibilities for 
the nature of this planet: (1) an ocean planet, (2) a bare rock 
with an iron content less than that of Earth, or (3) a terrestrial 
planet with a thin, low molecular weight atmosphere (e.g., 
H-He). The water world hypothesis is supported by a 
temperate equilibrium temperature for TOI-1452b of 298 + 
6K assuming an Earth-like Bond albedo (Ag =0.3), and 
between 226 and 326K for extreme Ag of 0.77 (Venus-like) 
and O (pure absorber). The insolation level of TOI-1452b is 
about 80% higher than the Earth (S = 1.8 + 0.2 S.), similar to 
Venus (S= 1.91 S,). 

Focusing on the first two possibilities, we used an MCMC 
approach (emcee; Foreman-Mackey et al. 2013) coupled to an 
interior structure model (Valencia et al. 2007) to obtain the 
distribution of mass fractions of iron and water that are 
consistent with the data. The details of this modeling can be 
found in Plotnykov & Valencia (2020). The HO layer is 
described by the equation of state from Hemley et al. (1987) 
and Stewart & Ahrens (2005). We report below the 16th, 50th, 
and 84th percentiles of the posterior distributions (available in 
Appendix C). 

From the chemical analysis of the star (Section 3.4), we 
obtain chemical ratios relevant to the planetary interior. 
Notably, TOI-1452 might have a slightly lower Fe/Mg weight 
ratio compared to the Sun (see Table 9), but in fair agreement 
with that of a sample of ~1000 M dwarfs from the APOGEE 
(DR16) catalog (Majewski et al. 2016; Ahumada et al. 2020) 
with known chemical ratios (see Figure 16). The APOGEE 
abundances are derived from high-resolution near-infrared 
spectroscopy (R ~ 22 000), with typical uncertainty on [Fe/H] 
and [Mg/H] of 0.02 dex. Obtaining the Si ratios for TOI-1452 
was difficult given the scarcity of spectral lines. In addition, we 
derive a C/O weight ratio consistent with the solar value. Our 
chemical abundance ratios for TOJ-1452 are summarized in 
Table 9. 

To infer the planet’s composition for scenarios | and 2, we 
can either make no assumptions on the refractory ratios, and 
thus, obtain all possible compositions that fit the mass—radius 
data; or, assume that the refractory ratios of planet and star are 
related, and use the star’s ratios as priors in the Bayesian 
analysis. Given that the refractory ratios of super-Earths seem 
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Figure 13. RV time series from SPIRou and IRD with the best-fit Keplerian + activity GP (orange curve), activity GP-only (red dashed curve), and Keplerian-only 
(blue curve) models overplotted. The residuals below show an overall agreement between the SPIRou errors and the RV dispersion, which is difficult to assess for IRD 


given the limited number of measurements. 


to span a larger range than that of stars (Plotnykov & 
Valencia 2020), we applied both methods here. In both cases, 
we kept the Mg/Si ratio in the planetary mantle the same as the 
star. This assumption should not affect the results considerably, 
given that the mantle minerals formed by different Mg/Si 
ratios have similar equations of state and thus, the relative 
content of Mg to Si is not constrained by planetary mass and 
radius data (Plotnykov & Valencia 2020). 

In the case where we make no assumptions, we obtain a 
water-mass fraction (WME) of 0.277072 and core-mass fraction 
(CMF) of 0.307279. These values translate to Fe/Si = 2.934 
and Fe/Mg = 54003 by weight. The planetary refractory 
ratios are particularly large primarily due to large uncertainty in 
planetary mass but also the degeneracy that ensues when 
considering water. The tail of the Fe/Mg distribution is long 
because this case allows for little to no mantle in the planet (see 
Figure C1). If instead, we assume the refractory ratios of the 
star as priors, we obtain lower water- and core-mass fractions: 
WME = 0.227013, CMF=0.18+0.06, resulting in Fe/ 
Si=1.340.4 and Fe/Mg=1.5 40.4. Both abundance sce- 
narios yield a nonzero, yet poorly constrained, WMF. The large 
uncertainty on both the CMF and WMF is rooted to the modest 
mass constraint (3.70). Better mass measurements are needed 
to confirm that TOI-1452b has a significant WMF. 

A different possibility (scenario 2 discussed above) is that this 
planet is a bare rock with no significant atmosphere, perhaps 
because it lost any acquired water through atmospheric evaporation 
during the high-insolation phase of the M dwarf host star (Barnes 
& Heller 2013; Bolmont et al. 2017). In this case, we constrain the 
physical model solutions but make no assumptions on the 
refractory ratios. The results show a CMF = 0.197018. with 


Fe/Si=1.07)2 and Fe/Mg=1.2')2 by weight. This CMF 
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indicates a planet not as dense as the Earth, with refractory ratios 
still consistent at the 1o level with those of its host star (see Table 9 
and Figure 16). However, the maximum a posteriori estimate of 
Fe/Mg (~0.8) for this bare rock scenario is lower than that of its 
host star. Forming planets that are iron poor with respect to their 
host star is difficult (Scora et al. 2020). Therefore, based on our 
current knowledge of planet formation and the chemical 
characteristics of the star, this scenario is less likely. 

A summary of this interior modeling is presented in Table 9 
and Figure 16. The posterior distributions for the three models 
(no assumptions, stellar priors and bare rock) are shown in 
Appendix C. 

Outside of these scenarios, our observations do not rule out 
other structures with a low molecular weight atmosphere 
(scenario 3) such as an Earth-like interior surrounded by an H— 
He envelope at T= 300 K containing ~0.1% of the total mass 
(up to 0.5% H-He at Io). One way to firmly break the 
degeneracy in planetary internal structures would be to 
characterize the atmosphere of TOI-1452b. 


5.2. Atmospheric Characterization Prospect 


TOI-1452b is a prime target for follow-up transit spectrosc- 
opy with JWST. The system is located near Webb’s continuous 
viewing zone (CVZ), more precisely at a few degrees (~10°) 
off the Northern CVZ, which means that it could be observed 
most of the year. Moreover, TOI-1452b is one of the few 
identified super-Earths in a temperate regime (TJ, between 200 
and 400K) orbiting a relatively bright star amenable to 
transmission spectroscopy observations (see Figure 17 and 
Table 10). The expected strength of the atmospheric signal is 
characterized by the transmission spectroscopy metric (TSM; 
Kempton et al. 2018), which is proportional to the host star’s J 
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Prior and Posterior Distributions of the Joint Transit RV Fit for Model M,.p;cp (Details in Section 4.3) Using Only the SPIRou Radial Velocities 


Parameter 


Prior* 


Posterior 


Description 


Fitted parameters 


px (gem) 

P (days) 

to (BJD - 2 457 000) 
ry 

"2 

K(ms_') 

e 

w (°) 

Agp (ms~') 
cp (days) 
‘spirou (m s_') 
OspiRou (mM s—') 
q1,TESS 

42,TESS 

Dress 

Mress 

Oress (ppm) 
41,PESTO 
42,PESTO 
Dpesto 
Mpesto 
Opesto (ppm) 


N(16.8, 1.97) 
U(11.0, 11.1) 
U(1691.4, 1691.6) 
uo, 1) 

Uo, 1) 
uo, 10) 

O (fixed) 

90 (fixed) 
LU0-?, 100) 
LUA0-?, 100) 
U(—33995, —33975) 
Uo, 10) 
uo, 1) 
uo, 1) 

1.0 (fixed) 
NO, 0.127) 
LUC, 10 000) 
Uo, 1) 
uo, 1) 
N(0.564, 0.05647) 
NO, 0.27) 
LUA, 10 000) 


16.8114 
11.06201 + 0.00002 
1691.5321 + 0.0015 

0.46 + 0.08 
0.0555 + 0.0014 
3.50 + 0.94 
0 
90 
4.5t79 
11.3443 
—33985 +2 
2.3413 
0.357073 
0.377933 
1.0 
-0.00032 + 0.00010 
15t7$ 

0.671051 

0.467) :38 


0.586 + 0.040 
0.00013 + 0.00018 


2287 + 67 


MuSCAT3 photometric parameters in Table 8 


Stellar density 

Orbital period 

Time of inferior conjunction 
Parameterization” for R,/R, and b 
Parameterization? for R,/R, and b 
RV semi-amplitude 

Orbital eccentricity 

Argument of periastron 
Amplitude of the GP 

Timescale of the GP 

SPIRou RV systemic component 
SPIRou RV extra white noise® 
TESS limb-darkening parameter® 
TESS limb-darkening parameter® 
TESS dilution factor 

TESS baseline flux 

TESS extra white noise 

PESTO limb-darkening parameter‘ 
PESTO limb-darkening parameter‘ 
PESTO dilution factor 

PESTO baseline flux 

PESTO extra white noise 


M, (Mz) 

p (gem) 

a (au) 

6 = (Rp/R.)” 

b 

i(°) 

Teg (K) 
[Ag = 0] 
[Ap = 0.3] 
[Ag = 0.77] 

S (Se) 


Derived parameters 

1.672 + 0.071 

4.82 + 1.30 
5.6118 

0.061 + 0.003 
3.09 + 0.16 
0.19 + 0.13 
89.77 + 0.16 


326 +7 
298 + 6 
226 + 5 
1.8+0.2 


Planetary radius 
Planetary mass 
Planetary bulk density 
Orbital semimajor axis 
Transit depth 

Transit impact parameter 
Orbital inclination 
Equilibrium temperature 


te 5 Tnsolation 


Notes. 


* U(a, b) is the uniform distribution between values a and b. LU (a, b) is the log-uniform (Jeffreys) distribution between values a and b. N(ju, a?) is the normal 


distribution with mean j: and variance 0”. 
> Parameterization from Espinoza (2018). 
© White noise term for single exposures within polarimetric sequences. 


{q1, 2} are linked to the quadratic limb-darkening coefficients {u;, u2} through the transformations outlined in Kipping (2013). 


Table 8 


Prior and Posterior Distributions of the MuSCAT3 Photometric Parameters for Model M,,p+cp (Details in Section 4.3) Using Only the SPIRou Radial Velocities 


Parameter 


Prior 


/ 


& 


rad 


of 


Fitted parameters 


o (ppm) 


uo, 1) 
Uo, 1) 
u(0.5, 1) 
N(O, 0.12) 
LUC, 10 000) 


0.601036 
y 
0.594032 
0.71 £0.12 


—0.00141 + 0.00036 


2748 + 319 


0.225.153 
0.377933 
0.887002 
—0.00106 + 0.00016 
463" 


Derived parameters 


Ouncorr. (ppt) 


+0.40 
2.18736 


9.727021 


L &s 
0.461035 0.241072 
yl 
0.244028 0.277033 
0.967).83 0.85 + 0.07 
—0.00130 + 0.00014 —0.00131 + 0.00013 
18371434 297327 
2.94 + 0.17 2.64 + 0.20 
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Figure 14. Phase-folded SPIRou RV curve, with systemic velocity and activity 
GP removed. Binned RV measurements (0.1 phase bin) are marked with black 
circles. The best-fit Keplerian models (blue) are depicted with a solid curve 
(circular orbit) and a dashed curve (eccentric orbit). The residuals of the 
circular fit are presented below. The circular solution yields smaller residuals 
rms (3.58m s_!) compared to the eccentric model (3.66 m s}); yet both 
Bayesian evidences are indistinguishable (see Figure 12). 
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Figure 15. Mass-radius diagram of exoplanets (gray points). Only planets with 
mass and radius known with precision better than 30% are shown. Various 
theoretical composition curves are plotted, using interior structure models from 
Valencia et al. (2007), Valencia et al. (2013), Guillot & Morel (1995). All 
models assume no atmosphere, except the H—He (red dashed lines) models that 
correspond to Earth-like interior surrounded by a solar composition envelope at 
T = 300 K. The mass and radius posteriors (green | and 2c contours) of TOI- 
1452b are consistent with either a water-rich interior with ~25% HO by mass, 
a pure rock (rocky threshold radius curve including relevant phase transitions), 
or a terrestrial planet with a thin hydrogen envelope of ~0.1% H—He by mass. 


magnitude and the planet’s atmospheric scale height. Figure 17 
displays the TSM as a function of the equilibrium temperature 
for known small exoplanets with available mass measurements. 
The sample is restricted to systems with well-determined 
masses (relative uncertainty <30%) as a constraint on surface 
gravity is essential to correctly interpret the transmission 
spectrum of an exoplanet (Batalha et al. 2019). The temperate 


17 


Cadieux et al. 


subset of Figure 17 is detailed in Table 10. The TSM of TOI- 
1452b (39.9) is similar to well-known temperate systems such 
as LHS 1140b (50.0) and K2-18b (40.8), while being 60% 
below the highest listed target in this subset, L231-32 d (TOI- 
270 d, 104.0). All seven host stars in Table 10 have Tere < 4000 
K (or average Tz; of 3225 K), confirming the high interest of M 
dwarfs for planetary atmospheric characterization. Note that the 
high-value target L 98-59 d (TSM above 200; Cloutier et al. 
2019; Demangeon et al. 2021) was just barely excluded from 
Table 10 due to its T.g (409 K) being slightly above our 400 K 
cut. Our subset also excludes intriguing planets with plausible 
temperate environment but deprived of mass measurements (or 
imprecise mass), such as TOI-700 c (TSM = 79.7; Gilbert et al. 
2020), TOI-1266c (TSM = 48.8; Demory et al. 2020) and K2- 
3c (TSM = 25.5; Damasso et al. 2018). The TOI-1452 system 
is a unique target to explore the atmospheric properties of 
temperate planets within the radius valley. This paper provides 
the first mass determination needed for the interpretation of 
future transmission spectra. 


5.3. Implications for the Emergence of the M Dwarf Radius 
Valley 


Planets on either side of the radius valley differ by their 
composition, typically “rocky” for the smaller super-Earths, 
and “gaseous” for the larger mini-Neptunes. This transition 
occurs as a consequence of a varying envelope mass fraction: 
adding an H—He envelope up to a few percents of the total mass 
of a planet essentially doubles its observable radius (Lopez & 
Fortney 2014; Chen & Rogers 2016). Thermally driven 
atmospheric escape processes such as photoevaporation (Owen 
& Wu 2013; Lopez & Fortney 2014; Owen & Wu 2017; Lopez 
& Rice 2018; Wu 2019) and core-powered mass loss (Ginzburg 
et al. 2018; Gupta & Schlichting 2019, 2020) have been 
proposed as radius valley emergence mechanisms. In these 
models, super-Earths and mini-Neptunes originate from the 
same population of planets that form with an extended H—He 
envelope around an Earth-like core, with the population of 
rocky super-Earths emerging after losing their primordial 
atmospheres to hydrodynamic escape. Another possible 
scenario is to assemble rocky super-Earths at late times after 
most or all of the gas has been dissipated from the 
protoplanetary disk (Lee et al. 2014; Lopez & Rice 2018; 
Lee & Connors 2021). The two classes of planets would form 
on different timescales, resulting in a bimodal distribution 
without relying on any subsequent atmospheric escape. 

Each of the aforementioned mechanisms predicts that the 
rocky-to-gaseous transition (Ryatey) varies with parameters 
such as the orbital period P and stellar mass M,. Photo- 
evaporation, core-powered mass-loss and gas-poor accretion 
models predict a negative slope in the R,—P space: 
Ryatiey x P-°?>~°!> (Qwen & Wu 2017; Lopez & Rice 2018; 
Mordasini 2020), Ryatiey xP °'' (Gupta & Schlichting 2019), 
and Rago (Lee & Connors 2021). Conversely, the 
formation of super-Earths strictly by the merging of planetary 
embryos in a gas-depleted environment, analogous to the 
formation of terrestrial planets in the solar system, would 
produce a positive slope (Ryatiey & pe", Lopez & Rice 2018). 
One way to test the proposed models is to compare these 
predictions to the real population of exoplanets. 

From occurrence rate calculations of small close-in planets 
around Sun-like stars, Martinez et al. 2019 measured a 
d log Ryattey /d log P = —0.11 + 0.02, consistent with 
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Table 9 
Chemical Ratios by Weight for the TOI-1452 System 
Chemical TOI-1452 Sun* TOI-1452b TOI-1452b TOI-1452b 
Ratios No assumptions Stellar priors Bare rock 
Fe/Mg 1.484036 1.83 + 0.25 3.4483 1.5104 1.255% 
Mg/Si 0.867037 1.06 + 0.13 0.86" 0.86" 0.86° 
C/O 0.487944 0.41 + 0.07 : ' 
CMF 0.307979 0.1810-06 0.197918 
2 
WMEF 0.2770-20 0.2249: . 
Notes. 


* Photospheric abundance ratios from Asplund et al. (2009). 
> The Mg/Si ratios of the planets are fixed to the star TOI-1452 ratio. 
© The interior models assume no carbon compounds. 
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Figure 16. Fe/Mg distribution for TOI-1452b according to different model 
assumptions (Table 9) compared to that of its host star (black). We used kernel 
density estimate to draw the probabilities from the posteriors. We also include 
the M dwarfs distribution of Fe/Mg from the APOGEE database (gray 
histogram, sample size of ~1000; Majewski et al. 2016; Ahumada et al. 2020) 
for comparison. 


thermally driven mass loss and gas-poor formation. Using a 
similar methodology but for planets around low-mass stars with 
Tere << 4700 K, Cloutier & Menou (2020) obtained a 
d log Ryatey /d log P = 0.058 + 0.022, where the positive 
sign suggests that the gas-depleted formation of super-Earths 
may be dominant around M dwarfs. These distinct slope 
measurements carve out regions in the R,—P parameter space 
where the models make opposing predictions regarding the 
bulk composition of a planet (i.e., either rocky or gaseous). 
This framework to test radius valley emergence models around 
M dwarfs was introduced in Cloutier & Menou (2020) and has 
since been applied to a number of transiting planets (TOIs 
776 b, Luque et al. 2021; 1235b, Cloutier et al. 2020b; 1634b, 
Cloutier et al. 2021; 1685b, Bluhm et al. 2021). 

Figure 18 presents the period—radius diagram for exoplanets 
around M-dwarf hosts (Tere< 4000 K). Each planet is color 
coded by its bulk density relative to the Earth-like structure 
model of Valencia et al. (2007; see Figure 15). TOI-1452b sits 
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on or slightly above the empirical valley of Cloutier & Menou 
(2020), while being considerably below the slope measured by 
Martinez et al. (2019), scaled down to match the median stellar 
mass of the Cloutier & Menou (2020) sample (using Equation 
(11) therein). The locus of TOI-1452b in Figure 18, combined 
with our density estimate, are incompatible with the photo- 
evaporation and core-powered mass-loss models. However, the 
likely intermediate nature of TOI-1452b cannot strongly 
support the alternative gas-depleted formation scenario either 
as the dominant mechanism for the emergence of the M-dwarf 
radius valley. A volatile-rich interior for TOI-1452b could 
indicate a different formation pathway, e.g., one without 
significant gas accretion during the disk lifetime. 

Figure 18 also highlights three other systems presenting 
similarities with TOI-1452b, namely, TOI-1235b (Cloutier et al. 
2020b), L 98-59 d (Cloutier et al. 2019; Demangeon et al. 2021), 
and LHS 1140b (Dittmann et al. 2017; Lillo-Box et al. 2020). All 
four planets have a similar size, while spanning a large interval in 
periods. TOI-1235b (P = 3.445 days, Rp = 1.738 Re) and LHS 
1140b (P = 24.737 days, Rp = 1.635 Ra) have densities compa- 
tible with bona fide super-Earths; their position in Figure 18 
indicates that they are probably examples of the largest terrestrial 
planets that can be assembled around M dwarfs without accreting a 
substantial hydrogen envelope. On the other hand, L 98-59 d 
(P=7.451 days, Rp = 1.521 Re) is a likely water-rich (~30%) 
planet that may be approaching, like TOI-1452b, the minimum 
size for volatile-rich objects around a low-mass star. These four 
systems constitute benchmarks for understanding the formation 
and evolution of planets within the radius valley. 


6. Summary & Conclusion 


This paper reports the discovery and characterization of the 
transiting temperate exoplanet TOI-1452b. A joint analysis of 
transit observations from TESS and other ground-based 
telescopes combined with radial velocity measurements from 
SPIRou and IRD, yields a mass of 4.82 + 1.30 Mz and a radius 
of 1.672 + 0.071 Rz. These physical parameters are consistent 
with either a rocky world with a Fe/Mg ratio similar to the host 
star (Fe/Mg =1.2*}2 by weight), a water-rich interior (either 
22+11% HO by weight, if stellar priors are assumed for the 
planetary refractory ratios, or 27*72% H,O if no assumptions 
are made) or a terrestrial planet surrounded by a <1% H-He 
atmosphere. Orbiting its M4 host star (Ter = 3185 +50 K) 
every 11.06201 +0.00002 days, the planet receives about 
twice as much radiation than the Earth (S=1.8+0.2 Sq), 
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Table 10 
TSM for a Subset of Well-characterized Small Exoplanets in a Temperate Equilibrium Temperature Regime (200 K < T., < 400 K) 
Planet P My Rp Teg [Az = 9] J Tete M, R, TSM References 
(days) (Mz) (Re) (K) (mag) (K) (Mo) (Ro) 

L231-32d 11.380 4.78 2.133 388 9.10 3506 0.39 0.38 104.0 (1) 
TOI-1231b 24.246 15.4 3.65 331 8.88 3553 0.48 0.48 97.6 (2) 
LTT 3780c 12.252 8.6 2.30 353 9.01 3331 0.40 0.37 72.6 (3) 
LHS 1140b 24.739 6.38 1.635 214 9.61 2988 0.19 0.21 50.0 (4) 
K2-18b 32.940 8.63 2.610 279 9.76 3457 0.50 0.44 40.8 (5) 
TOI-1452b 11.062 4.82 1.672 326 10.60 3185 0.25 0.28 39.9 (6) 
TRAPPIST-1 b 1.511 1.374 1.116 399 11.35 2566 0.09 0.12 28.6 (7) 
TRAPPIST-1 d 4.049 0.388 0.788 287 11.35 2566 0.09 0.12 25.6 (7) 
TRAPPIST-1 c 2.422 1.308 1.097 341 11.35 2566 0.09 0.12 24.3 (7) 
TRAPPIST-1 e 6.101 0.692 0.920 251 11.35 2566 0.09 0.12 19.9 (7) 
LHS 1140c 3.777 1.76 1.169 400 9.61 2988 0.19 0.21 18.7 (4) 
TRAPPIST-1 f 9.208 1.039 1.045 218 11.35 2566 0.09 0.12 17.0 (7) 
Note. 


References: (1) Van Eylen et al. (2021). (2) Burt et al. (2021). (3) Cloutier et al. (2020a). (4) Lillo-Box et al. (2020). (5) Benneke et al. (2019). (6) This work. (7) Agol 


et al. (2021). 
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Figure 17. TSM (Kempton et al. 2018) as a function of the planetary 
equilibrium temperature for small exoplanets (Rp < 4 Rg) with well-established 
masses (error below 30%). The green region corresponds to an arbitrary 
temperate temperature interval (200 K < T., < 400 K). TOI-1452b is among 
the best temperate targets for follow-up transit spectroscopy with JWST. 


corresponding to a blackbody temperature of 326 +7 K. The 
results of our interior modeling and the fact that the planet 
receives modest irradiation make TOI-1452b a good candidate 
water world. 

TOI-1452b is a prime target for upcoming atmospheric 
characterization efforts with JWST, featuring a high transmis- 
sion spectroscopy metric compared to other known temperate 
exoplanets. Transit spectroscopy observations with JWST 
should reveal the true nature of this intriguing exoplanet lying 
within the radius valley, whether this is a rocky world or one 
with a volatile envelope. Being observable with JWST most of 
the year, TOI-1452b is a unique system for studying exoplanets 
at the transition between super-Earths and mini-Neptunes. 


We acknowledge the use of public TESS Alert data from 
pipelines at the TESS Science Office and at the TESS Science 
Processing Operations Center. This research has made use of 
the Exoplanet Follow-up Observation Program website, which 
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Figure 18. Period—radius diagram of exoplanets around low-mass stars 
(Tepe < 4000 K) with mass and radius known with precision better than 30%. 
The TRAPPIST-1 (T-1) planets are represented by triangles. Empirical valleys 
for low-mass stars from Cloutier & Menou 2020 (CM20), consistent with gas- 
depleted formation, and from Martinez et al. 2019 (M19), compatible with 
thermally driven atmospheric mass loss and scaled down to represent the stellar 
mass population of CM20. TOI-1452b lies below M19 and has a small bulk 
density compared to a terrestrial planet of the same mass (~/Pgarth—tike © 0.76). 
This is inconsistent with photoevaporation and core-powered mass-loss 
predictions. 
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Appendix A 
TESS Light Curve 


We present the TESS multisector PDCSAP light curve in 
Figure Al, with the exception of sectors 14 and 21, which were 
previously shown in Figure 2. 
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Figure Al. Normalized PDCSAP light curve of TOI-1452 from sectors 15-26 (except 18), 40-41, and 47. Sectors 14 and 21 were previously shown in Figure 2. The 
blue data points highlight the epochs of the detected transits, while the red data points are outliers /stellar flares rejected by sigma clipping (3.5 o clip) or manually (in 
sectors 40, 41, or 47). The gaps in the light curve coincide with data downlink when TESS is close to perigee. The full light curve was modeled with a quasiperiodic 
Gaussian Process (green curve; details in Section 4.2). 
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Appendix B 
Supplementary Material Regarding the Joint Transit 
RV Fit 


In this appendix, we summarize the RV component of the 
joint transit RV models introduced in Section 4.3. The main 
RV parameters of the joint fits are reported in Table B1. 
All models and data sets detect TOI-1452b with a coherent 
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semi-amplitude K. Models with an activity GP (Mjop+cp and 
Miep+cp) produced the highest Bayesian log-evidence 
(Figure 12) and needed the least extra white noise (Cspyrou; 
O1rp). We ultimately adopted the results of the SPIRou only 
joint fit because adding the seven IRD measurements yields 
similar Bayesian log-evidences and requires a Oygp about three 
times the level of the planetary signal (see Table B1). 


Table B1 
RV Component of the Joint Transit RV Fit for Different Models (M) and Data Sets 
Parameter lcp lcp+GP lep lep+GP 
SPIRou only 
K(ms7!) 4.2+0.9 3.540.9 4.7 £0.9 3.6 + 0.9 
e = 0.20 + 0.09 0.124012 
Agp (m s_') 4.5ttS 4.4tt9 
cp (days) 11.3439 ve 11.6724 
OspiRou (ms!) 49+0.9 23413 4.7108 22413 
SPIRou + IRD 
K(m s!) 41+0.9 3.52 0.9 4.7+0.9 3.6+0.9 
e a 0.19 + 0.09 0.124022 
Agp (m s~!) 473 4.6175 
lcp (days) 111th? wi 11.744}? 
Ospirou (ms!) 4.9+0.9 2.3413 4.7106 23805 
o1rp (m s~!) 13.6134 11.5434 14.2433 11.9737 


Notes. Mp: single circular orbit planet. Mop +Gp: 
single eccentric orbit planet and activity GP. 
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Appendix C 


Summary of Interior Parameters ae : : 
the interior, respectively, when no assumptions on _ the 


We present in this appendix the corner plots from the refractory ratios (Fe/Si and Mg/Si) are made and when the 
MCMC modeling of the planetary interior. Figures Cl and C2 planet follows the refractory ratios of its host star. Figure C3 
summarize the cases when water is allowed to be significant in presents the result of the bare rock scenario. 
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Figure Cl. No assumptions, corner plot summary for TOI-1452b interior parameters where the core-mass fraction (CMF) and water-mass fraction (WMF) are 
simulated quantities for a given planetary mass and radius. Chemical ratios Fe/Si and Mg/Si are derived quantities. The vertical red lines in the Fe/Si, Fe/Mg space 
are the mean stellar refractory ratios, and the vertical dashed lines represent 16th and 84th percentiles. 
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Figure C2. Stellar prior, corner plot summary for TOI-1452b interior parameters where we assume that the planet follows stellar refractory ratios for a given planetary 
mass and radius. Thus, the log-probability function is modified to include the restriction posed by the Fe/Mg ratio of the star and can be written as 
Fe/Mg ~ N(Fe./Mg,, oFe,/Mg,)- The red and dashed vertical lines represent the stellar ratios and 16th and 84th percentiles of the posterior, respectively. 
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Figure C3. Bare rock, corner plot summary for TOI-1452b interior parameters where we assume only rocky composition is possible (no water) for a given planetary 
mass and radius. Chemical ratios Fe/Si and Mg/Si are derived quantities. The red and dashed vertical lines represent the stellar ratios and 16th and 84th percentiles of 
the posterior, respectively. 
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Appendix D 
Radial Velocity Measurements 


We present the radial velocity measurements of TOI-1452 
from SPIRou and IRD in the online Table D1. 


Table D1 
SPIRou and IRD RV Measurements 


Instrument BJD - 2 400 000 RV (ms_') ory (ms!) 
SPIRou 59004.995291 —33983.32 8.39 
SPIRou 59005.006067 —33975.16 8.12 
SPIRou 59005.016836 —33980.21 8.25 
SPIRou 59005.027551 —33995.65 8.28 
SPIRou 59009.008115 —33975.18 8.02 
SPIRou 59009.018950 —33974.66 8.20 
SPIRou 59009.029725 —33975.25 7.96 
IRD 59118.766592 15.17 4.08 
IRD 59122.769496 17.14 3.86 
IRD 59122.787278 14.54 3.77 
IRD 59156.837456 —14.4 4.86 
IRD 59372.001844 —9.95 3.71 
IRD 59373.904952 2.00 5.70 
IRD 59390.850628 12.31 4.03 


Note. Table D1 is published in its entirety in the machine-readable format. A 
portion is shown here for guidance regarding its form and content. 


(This table is available in its entirety in machine-readable form.) 
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